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Abstract 

We re-fit for the neutrino mass-squared difference Am? in the threefold maxi- 
mal (ie. tri-maximal) mixing scenario using recent CHOOZ and SUPER-K data, 
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H taking account of matter effects in the Earth. While matter effects have little 

influence on reactor experiments and proposed long-baseline accelerator exper- 
iments with L ^ 1000 km, they are highly significant for atmospheric experi- 
ments, suppressing naturally I'e mixing and enhancing i/^ — Vj- mixing, so as to 
effectively remove the experimental distinction between threefold maximal and 
twofold maximal — Vr mixing. Threefold maximal mixing is fully consistent 
with the CHOOZ and SUPER-K data and the best-fit value for the neutrino 
mass-squared difference is Am? ~ (0.98 ± ) x 10^^ eV^. 
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1. Introduction 



New limits on neutrino oscillations from the CHOOZ reactor experiment []I| defini- 
tively rule out threefold maximal lepton mixing with Am^ ~ 10^^ eV^ [0. In the 
meantime however, the second-generation underground water-Cherenkov experiment 
SUPER-KAMIOKANDE has superseded the older KAMIOKANDE experiment § 
as regards measurements of the atmospheric neutrino anomaly, with the latest exper- 
imental fits to z/^ — Vr mixing [^] clearly favouring a much smaller value of Am^ than 
before 0. The aim of this paper is to underline the fact |0 that, in the light of such 
developments, threefold maximal mixing is per se far from excluded, in particular we 
emphasise, when terrestrial matter effects are taken into account. 

In this paper we re-fit for the neutrino mass-squared difference Am^ in the threefold 
maximal mixing scenario using the CHOOZ and SUPER-K data, taking full account 



of terrestrial matter effects which were hitherto neglected ^ |T0[ . This simplest of 
all possible mixing schemes then continues to explain the great majority of oscillation 
data just as previously claimed, but now with Am^ ~ 10~^ eV^. 



2. The Vacuum Scenario 



In vacuum in the threefold maximal mixing scenario, the mixing matrix relating the 
neutrino flavour eigenstates z/g, t'^, Ut to the neutrino mass eigenstates z/i, z/2, (with 
masses mi < m2 < respectively) takes the symmetric or 'democratic' form: 
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where here and throughout this paper we give directly the squares of the moduli 
of the mixing elements, in place of the mixing elements themselves. The mixing 
phenomenology is then completely determined by the two independent vacuum mass- 
squared differences Am^ > Am'^ (eg. Am? = m\ — Am'^ = m\ — m\). From the 
atmospheric data we have Am^ ~ 10~^ eV^ (see below) while from the solar data we 
found Am'^ < 0.9 x 10^^^ eV^ 0, so that in threefold maximal mixing, two masses 
are effectively degenerate and the spectrum of mass-squared differences is hierarchical. 

Neutrino oscillations violate lepton flavour conservation and as a function of prop- 
agation length L the matrix of normalised transition amplitudes Ain from a charged 
lepton state I to a charged lepton state /' is given directly by exponentiating the neu- 
trino mass matrix squared mm^ (= m?) in the flavour basis: A = exp{—im'^L/2E), 
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where E is the energy of the neutrino, assumed to be relativistic. The mixing matrix 
U comprising the normahsed eigenvectors of diagonahses the mass matrix thus: 
Wm'^U = diag ( ml, m\, m\ ). In terms of the mixing matrix the above matrix ex- 
ponentiation is readily achieved by first exponentiating the diagonal matrix of phases 
and then rotating back to the fiavour basis ('un-diagonalising') as follows: A = U diag 
( e**^^, e*"^^, e*"^^ )U\ where the phases (pi = mfL/2E for i = 1 ... 3. Thus any transi- 
tion amplitude Aiq may be written as the sum of three sub-amplitudes A^q = X^^i e'^'^' , 
where Xi^i = UinUi^, corresponding to the independent propagation of each of three 
neutrino mass eigenstates as usual. (Note that hypothetical 'sterile' neutrinos []Tl| 
play no role in the threefold maximal mixing scenario and are not considered here.) 

Survival and appearance probabilities P{1 I) = {Aul"^ and P{1 I') = are 
then given by: P{1 ^ I) = |A« + A^'> + A^^\' and P{1 I') = + Af} + Af)\\ 
respectively. In threefold maximal mixing the three sub-amplitudes corresponding to 
the three mass eigenstates are always of equal modulus 1/3. All vacuum survival 
probabilities P(/ — > /) and appearance probabilities P{1 I'), P{1' — > I) are indepen- 
dendent of fiavour (/, V = e, fi, r). In the limit Am'^ only one relative phase need 
be retained (eg. 03 — 02; where 02 = 0i for Am'^ = 0) and we are led to: 

P(l ^ /) = 5/9 + 4/9 cos{Am^L/2E) (2) 

and: 

P{1 I') = p{l' ^ /) = 4/9 - 4/9 cos(Am=^L/2E). (3) 

Averaging over a range of E (and/or L) such that oscillating terms no longer con- 
tribute, still in the limit Am'^ with two sub-amplitudes adding 'coherently', one 
has for general mixing: < P{1 I) > = {xj;P + Xi^Y + (x/f-*)^ where the x//-* are 
real (the matrices are hermitian), and: < P{1 ^ I') > = |x/,|^-hX/,^^|2-h|x/,^^|2 = 
2|x/,;''p, where the last equality relies on unitarity (X"*^ + X^ + X'^ = /, the identity 
matrix). In the threefold scenario, one expects then a first 'step' or 'threshold' at 
L ~ (Am^/2i?)^^, marking the descent from P{1 —»/) = ! to: 

< P(/ ^ /) > = (1/3 + 1/3)2 + (1/3)2 = 5/9 (4) 

and the rise from P{1 I') = P{1' ^ I) = to: 

< P{1 /') > = < Pi^' ^ > = 2 X (l/3)(l/3) = 2/9. (5) 

If Am'2 is non-zero there will be a second threshold at L ~ {Am''^ /2E)~^ beyond 
which all three sub-amplitudes add 'incoherently' so that: < P{1 —>■ I) > = (XiP)'^ + 
{Xlp f + (X;fV with: <P{l^l')>= |X;,;'|2 + + |X/,?|2, and in threefold 

maximal mixing: <P{l^l)> = < P{1 ^ /') > = < P{1' ^ /) > = 1/3. 
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For solar neutrinos, in the threefold maximal mixing scenario, matter effects in the 
Sun are expected to be rather unimportant []T2[ and the vacuum prediction Eq. (4) 



may be compared directly with the measured suppressions, at least for the gallium 



experiments. The same comparison for the SUPER-K and HOMESTAKE \H 
solar results, (correcting for the neutral current contribution to the ue ue rate in 
SUPER-K) is limited by de facto uncertainties in the ®B flux 

In the atmospheric experiments the initial beam comprises both z/^ and Ue, and 
account must be taken of z/^ ^ Ue transitions, which tend to compensate the survival 
rates in Eq. (4). In the approximation that the flux ratio at production 0(z/^)/0(z/e) = 
2/1 the effective z/^ suppression in vacuum becomes: 

5/9 + 1/2 X 2/9 = 2/3 (6) 

while the effective suppression becomes: 

5/9 + 2 X 2/9 = 1 (7) 

so that the z/g rate is perfectly compensated in this case. The initial flux ratio is indeed 
0(z/^)/0(z/e) ~ 2/1 for ^ 1 GeV (increasing with E thereafter). 

While this concludes our discussion of the vacuum scenario, as we turn to consider 
the influence of matter effects it will prove useful to denote vacuum quantities with 
argument (0) for vacuum explicitly, thus we define: m^(0) = m^, U{0) = U, z/i(0) = Ui, 
mi(0) = nii, Am2(0) = Am^, Am'2(0) = Am'^ etc. 



3. Terrestrial Matter Effects 



In any scheme with z/g mixing (eg. Eq. 1), matter effects in the Earth, in particular 
for atmospheric neutrinos, can be very significant indeed. Matter effects result from 
the forward scattering of z/g from electrons in the matter, modifying the mass matrix 
by adding a term to the (e, e) entry of the vacuum mass matrix, proportional to the 
matter density p. Basic trends in matter phenomenology can usefully be anticipated 
by noting that in the limit of infinite density {p oo), the becomes the heaviest 
neutrino, which being then itself a mass eigenstate, completely decouples in the mixing. 
In general clearly, before this limit is approached, masses, mixings etc. are functions of 
the matter density, hence we have: m^(p), U{p), Vi^p), mi{p), Am^(p), Am'^(p) etc., 
or more directly functions of the number density Ng of electrons in the matter. 

In the present analysis matter effects are incorporated using the general 3x3 
numerical program described in Ref. [Q, where for a hierarchical spectrum Am^(O) ^ 
Am'^(O), the formulae of Ref. [|l^ apply- Our numerical calculation proceeds in steps 
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AL = 100 km through the Earth, with the density and composition of the Earth input 
as a function of depth from a recent tabulation |]TB[. Before giving our full numerical 
results in any detail however, in this section we sketch briefly the main trends predicted 
analytically in threefold maximal mixing with Am'2(0) < Am^{0), function of 
increasing matter density, or equivalently as a function of increasing E. 

We identify three significantly distinct density regimes where the mixing phe- 
nomenolgy remains essentially constant over a wide range of scales. 

Firstly, clearly at very low matter densities \/2GNe <^ Am'^(0)/-E (where G is the 
Fermi constant) vacuum mixing remains essentially valid. Conceivably if Am'^(O) is 
extremely small or zero (Am'^ ^ 10~^° eV^) this regime may not be realised in nature, 
even for neutrinos from distant supernovae propagating in deep space {N^ ~ 10^ m~^). 

Next, at 'intermediate' densities ie. for Am''^{0)/E < V2GNe < Am'^{0)/E, mat- 
ter effects lift the effective degeneracy between the two light neutrinos, such that in 
matter for neutrinos: z/i(p) — ^^2(0)] /a/2 and z/2(p) WiiO) + 1^2(0)] /V^ (up 

to phases), where the lighter matter mass eigenstate i^iip) has zero z/g content. The 
mixing matrix is thereby deformed in matter as follows: 
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with the smaller mass-squared difference ATn'^{p)/2E ^ 2/3 x \/2GNe in this case. 
The factor 2/3 multiplying the matter mass scale reflects the fact that only two out of 
the three relevant 2x2 sub-determinants of the vacuum mass matrix are modified by 
the matter term. (Symmetrically, this same factor of 2/3 will appear again multiplying 
the vacuum mass scale in the high density limit, see below). In the Earth's mantle 
(p ~ 5 g cm-3) the matter mass-scale ^ is v^GiVe ~ 0.38 x 10"^ eVVGeV, and 
the corresponding length-scale is {\/2GNe)^^ ~ 1040 km/rdn. Thus the limit Eq. (8) 
should be physically relevant for atmospheric neutrinos, assuming Am^(O) ~ 10^'^ 
eV^, for neutrino energies E 1 GeV. 

The matrix on the right-hand-side of Eq. (8) is one of the matrices with the 
maximally mixed introduced in Ref. ||^. For L ^ {\/2GN(,)~^ the vacuum mixing 
predictions Eqs. (2)-(7) are reproduced. For L > {■\/2GNe)~^ ^ ie. beyond the 'matter 
threshold', the i/^ survival probability averaged over E and L, is from row 2 of Eq. (8): 

< P(/i -> /i) > = (1/2)2 ^ (i/g)2 ^ (^/3)2 ^ 
while the average z^e survival probability is (from row 1): 

< P(e ^ e) > = (0)2 + (2/3)2 + (1/3)^ = 5/9 (10) 
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just as it was after the first threshold, defined by Am^, in the vacuum analysis above. 
Likewise the corresponding i/^ v^, appearance probabilities are: 

< e) > = < P(e ^ ^) > = (0)(l/2) + (2/3)(l/6) + (l/3)(l/3) = 2/9 (11) 

(from rows 1 and 2) again as in the vacuum analysis. 

For atmospheric neutrinos, in the approximation that the initial flux ratio (^ivp) / 
(pi^e) = 2/1, the effective suppression becomes: 

7/18 + 1/2 X 2/9 = 1/2 (12) 



replacing the vacuum prediction of 2/3 from Eq. (6) 0§[|ig], for L > {^/2GNe)-^. 
The Ve rate however remains perfectly compensated just as it was for / (f){Vf.) = 
2/1 in the vacuum analysis and Eq. (7) applies. Note that these final suppression 
factors are exactly the same as would be expected in twofold maximal z/^ — z/^ mixing. 

Finally, for sufficiently high matter densities Am'^(0)/E <^ \/2GN(. (or equivalently 
for sufficiently high neutrino energies), threefold maximal mixing tends physically to 
twofold maximal — mixing with Am^(O) fixing the smaller mass-squared differ- 
ence Am'^(oo) ^ 2/3 X Am^(O) in matter, and with the larger matter mass-squared 
difference fixed by the matter mass- scale Am2(oo)/2E ^GN^: 



y\{p) y2{p) y?,{p) i^i(oo) z/2(oo) 2/3(00) 



e 



I . 2/3 1/3 \ e 
1/2 1/6 1/3 1 — ^ /i 
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/ . .1 

1/2 1/2 . I . (13) 
V 1/2 1/2 



As foreseen, with the matter interaction affecting z/g directly, v^, —>■ z/3 asymptotically. 
Thus the higher frequency oscillations have zero amplitude whereby, for example, the 
z/t- appearance probability in a high energy z/^ beam becomes simply: 

P(^ ^ r) = 1/2 - l/2cos[2/3 x Am^{0)L/2E]. (14) 

For L — i> vacuum predictions must always be reproduced [§. Clearly Eq. (14) 
coincides with Eq. (3) for L — > 0. For L — > 00 we have < P(/i r) > ^ 1/2, so 
that appearance may be said to be enhanced by matter effects, cf. Eq. (5), and 
correspondingly < P{fi -h^ yu) > — > 1/2 in that case. As the z/g decouples completely 
P{e — >■ e) — i> 1 with P(/x — e) = P{e — /i) ^ in the limit, so that in the 
atmospheric experiments there are no significant compensation effects expected at the 
highest energies, despite the fact that the production flux ratio 0(z/^)/0(z/e) is becoming 
large, 0(z/^)/0(z/e) > 3/1 for ^ > 10 GeV. 

The important conclusion is that threefold maximal mixing with terrestrial mat- 
ter effects exactly mimics twofold z/^ — z/,- mixing for the atmospheric neutrino rates 
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both at the lowest energies and at the highest energies, but for somewhat different 
reasons in the two cases, as we have seen. Observable differences between threefold 
and twofold mixing are expected at best only in a window of 'intermediate' energies, 
where (piv^) / (pive) > 2/1 but before the limit Eq. (13) is approached. Of course our 
numerical results reported below incorporate the full expected energy and zenith-angle 
dependence of the flux ratio at production 0(z/^)/0(z/e) |TB[, with appropriate averag- 
ing over detailed u/u differences, neglected in the above discussion. 

4. Multi-GeV Zenith-Angle Distributions 

When oscillations are not individually resolved the neutrino mass-squared difference 
is determined by the location of the corresponding threshold on the L/E scale (the 
location of the 'matter threshold' is of course energy- independent, and is arguably 
better studied as a function of L rather than L/E). For atmospheric neutrinos the 
neutrino flight-path L is related to the zenith angle 6 of the neutrino by 

L = ^Rl cos2 e + 2R^H + m-R^ cos 6 (15) 

where Rq, ~ 6380 km is the radius of the Earth and ~ 20 km is the effective height 
of the atmosphere. In the water-Cherenkov experiments the neutrino zenith angle 
is estimated using the measurement of the outgoing charged-lepton, which will best 
correlate with the neutrino direction at higher energies, so that we expect the most 
reliable information to come from the so-called multi-GeV sample (ie. events with 
charged lepton momentum p ^ 1.3 GeV for electrons and p ^ 1.4 GeV for muons). 

Fig. 1 shows the measured zenith angle distributions for multi-GeV events, com- 
bining the SUPER-K ^ and KAMIOKA data for maximum statististical weight. 
There is clear evidence for a step (or 'threshold') as a function of cos^ with an ap- 
proximate 50% suppression of the /i-like events for cosO < —0.2. At the same time 
no such effect is apparent in the corresponding distribution for e-like events. Fig. lb. 
Since a full monte-carlo simulation including detector effects etc. is beyond the scope 
of the present paper, we have wherever possible in this analysis made appropriate use 
of the expected event rates for no oscillations given by the experimenters themselves. 
Thus in Fig. 1 we plot the ratio of observed to expected events in preference to the 
event rates themselves, in order to minimise dependence on detector acceptance etc. 
The plotted ratio in each case is normalised to the threefold maximal prediction with 
matter effects for Am? = 0.98 x 10"'^ eV^ (which is the overall best-fit Am^ in threefold 
maximal mixing, see Section 5 below) so that it is only the shape of the zenith-angle 
dependence which is being tested here, with no reliance on absolute fluxes. 
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In Fig. 1 the solid curves represent the threefold maximal mixing predictions com- 
puted taking full account of terrestrial matter effects and compensation effects for 
Am^ = 0.98 X 10~3 eV^ with the dotted curves showing the corresponding vacuum 
predictions. The dashed curves represent twofold maximal z/^ — Vr mixing for the 
same value of Arn^. All the above curves are averaged over the multi-GeV energy 
distribution and incorporate angular smearing. We take gaussian smearing in angle 
around the neutrino direction with a width fixed from existing v/v data |]21]], where the 
width falls with energy proportional to 1/ \/E. The mean angle between the neutrino 
direction and the charged-lepton direction is then 21° for the multi-GeV sample (and 
43° for subGeV events with p > 400 MeV/c, see Section 5 below). For comparison, 
angular smearing resulting from the Cherenkov method itself has the same energy de- 



pendence, but is an order of magnitude smaller |19] and may be neglected. Particle 
mis-identification effects, ie. mis-classification of yU-events as e-events and vice versa. 



believed to become important only at the percent level |2^, are also neglected in our 
analysis. As is clear from Fig. 1 (and see also Table 1, Section 5 below) both threefold 
maximal mixing with matter effects and twofold maximal z/^ — Ur mixing are consistent 
with the measured multi-GeV zenith-angle distributions. 

The predicted up/down ratios (U/D)^ and {U/D)f, for multi-GeV /x-like and e-like 
events respectively, are shown in Fig. 2, as functions of Aim? , for each of the various 
mixing scenarios above. The up/down ratios are defined from Fig. 1 as the ratio of 
up to down rates with | cos 6*1 > 0.2 and are expected ||2^ to be relatively free of 
systematic effects, eg. flux uncertainties. At multi-GeV energies the initial flux ratio 
0(z/^)/0(z/e) > 2/1, and in the case of threefold maximal mixing one expects to see 
o?;er-compensation of the Ve rate, ie. {U/D)e > 1, cf. Eq. (7), and under-compensation 
of the rate, ie. {U/D)^ < 1/2, cf. Eq. (12) (with smearing {U/D)^ ~ 1/2, Fig. 2), 
at least for Am^ ^ 10^'^ eV^. As is clear from Fig. 2, for Am? ^ 10^'^ eV^ over- 
compensation of the Ue rate is suppressed relative to the vacuum prediction as a result 
of terrestrial matter effects, with {U/D)e approaching unity as Am^ 0, simulating 
z/^ — Ur mixing as discussed above. 

In Fig. 2 the data points with error bars represent the measured up/down ratios 
based on the combined SUPER-K and KAMIOKA data: (U/D)^ = 0.53± and 
{U/D)e = 0.99± J-^y, these being plotted arbitrarily at Am^ ~ 0.98 x lO'^ eV^ and ex- 
tended by the shaded bands. The up/ down ratio for all events {U / -D)e+^ = 0.68±0.05 
(not shown) is less incisive, but has the advantage of being independent of particle 
misidentifiaction effects. Clearly the data on the up/down ratios are consistent with 
either twofold maximal — Ur mixing, or with threefold maximal mixing with inclu- 
sion of terrestrial matter effects, for Am? ~ 10"'^ eV^. 
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5. The Sub-GeV Data, R-values and Upward Muons 



The KAMIOKA/SUPER-K collaborations have also given zenith angle distributions 
for sub-GeV events (lepton momentum p < 1.3 GeV), as well as values for the double 
ratio R = {ft/e)DATA/ {fJ'/€.)MC for sub-GeV and multi-GeV events, and independent 
data on upward muons. For some of these data-sets there are significant systematic 
effects to be considered, against the gain in statistiscal weight. 

Fig. 3 shows the observed zenith angle dependence for the higher energy {p > 400 
MeV) subset of the sub-GeV sample in the SUPER-K experiment. The dashed curve 
is the prediction of threefold maximal mixing for Am^ ~ 0.98 x 10~^ eV^ including 
terrestrial matter effects and compensation effects, showing clearly the 'matter thresh- 
old' and the matter-induced oscillations, but neglecting angular smearing. The solid 
curve includes the expected effect of angular smearing as discussed previously, and 
which is now very significant. A fit to the data of Fig. 3 in threefold or twofold mixing 
yields a best-fit Am^ ~ 2 x lO"^ eV^ (xVDOF ~ 5.0/7, CL~ 66%). The zenith-angle 
data for the lower energy {p < 400 MeV) subset of the sub-GeV data are not included 
in this analysis. Angular smearing and geomagnetic effects are expected to domi- 
nate the zenith-angle dependence at such low energies, where very little zenith-angle 
dependence is seen and essentially no Am^ information survives. 

The integrated ratio of ratios R = {fi/ e) data/ ifJ'/ mc is independent of smearing 
effects, but dependent on flux uncertainties which affect the flavour ratio 0(z/^)/0(z/e), 
particularly since e-like and /i-like events have different Cherenkov thresholds. The 
latest i?-values given by the SUPER-K experiment are: R = 0.67 ± 0.02 ± 0.05 for 
the sub-GeV sample and R = 0.66 ± 0.04 ± 0.08 for the multi-GeV sample, where 
the first error is statistical and the second is the systematic error in each case, as 
given by SUPER-K. These results are shown in Fig. 4 by the data points and shaded 
bands, together with the predictions from threefold maximal mixing (solid curves) and 
twofold z/^ — Ur mixing (dashed curves), plotted as a function of Am^. Clearly within 
the combined statistical and systematic errors the consistency with Am^ ~ 10"'^ eV^ 
is perfectly acceptable. It is perhaps worth commenting here that threefold maximal 
mixing goes some way to resolving the problem of the near-equality of the multi-GeV 
and sub-GeV i?- values raised by LoSecco |^ . 

In the underground detectors upward muons result from neutrino interactions in 
the rock surrounding the detector |2^ [^. Fig. 5 shows the double ratio S of stopping 
to through muons for data divided by monte-carlo: S = 0.59 ±0.06 ±0.08 as measured 
by SUPER-K |^ (data point and shaded band). The solid curve is our threefold 
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maximal mixing prediction with terrestrial matter effects included, plotted as a func- 
tion of Am^, while the dashed curve corresponds to twofold maximal z/^ — Vj. mixing. 
The shift of a factor of 2/3 in Am? between the two curves as expected at high ener- 
gies, see Eq. (14), is apparent, with threefold maximal mixing showing a slightly lower 
minimum value for S", due to terrestrial matter effects. The relatively low measured 
value of S seems to point here to Am^ ~ 10^^ eV^ in either threefold or twofold mixing. 

6. The CHOOZ Data and Overall Fit 

The CHOOZ experiment |]I| measures the survival probability P(e — e) for z7g by 
comparing the observed to expected rates for the reaction z/g + p e+ + n at a 
distance L ~ 1 km from the CHOOZ reactor site. The initial result for P(e e) 
is shown in Fig. 6, plotted as a function of antineutrino energy related to the 
measured positron energy E^hy E = 1.8 + i?e MeV. The data so far are consistent 
with P(e — e) = 1, and no oscillation signal is claimed. 

At such short pathlengths L <^ (\/2GiVe)~^, matter effects can be safely neglected, 
and in Fig. 6 the solid curve is the threefold maximal mixing prediction calculated 
from the vacuum formula Eq. (2) for Am? = 0.98 x 10^'^ eV^. Twofold maximal 
^'/i ~ mixing predicts P(e — e) = 1 as shown by the dashed line, independent of 
matter effects and independent of E (and Am^). The best-fit to the data of Fig. 6 in 
threefold maximal mixing is: Am^ = 0.60 x lO^^ eV^ (xVDOF = 10.6/9, CL= 30%), 
although clearly the data are also fully consistent with Am^ 0, (x^/DOF 11.0/9, 
CL= 28%), or equivalently with — mixing. 

Our overall fit for the neutrino mass-squared difference Am? is based on the atmo- 
spheric data already discussed, taken together with the CHOOZ data for P(e ^ e) as a 
function of E. Specifically we compute the total summed over all the bins displayed 
in Figs. 1, 3, 4, 5, 6 as a function of Am^ for twofold maximal z/^ — Vr mixing and for 
threefold maximal mixing with matter effects. There are 33 data points in total, but 
4 normalisation factors to be determined in each case (one for each zenith-angle plot) 
in addition to Am^, so that there are 33 — 5 = 28 DOF in each fit. 

The results of the overall fits are shown in Fig. 7. For twofold maximal — 
mixing (dashed curve) there is a broad minimum in extending over the range 
Am^ ~ 7 X 10"^ - 4 X 10"^ eV^ with the best-fit Am^ ~ 2.2 X 10 ^ eV^ corresponding 
to an excellent fit (x^/DOF = 18.7/28, CL = 91%). For threefold maximal mixing 
(solid curve) the mass-squared difference is highly constrained by the CHOOZ data 
giving a much narrower minimum in x^, and a significantly more restricted Am?- 
range: Am? ~ (0.98 ± q'23) ^ 10^'^ eV^, however still corresponding to a very good 
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fit (x^/DOF = 25.4/28, CL = 61%). A full breakdown of the x^-contributions coming 
from the various data-sets is given in Table 1. 

While both threefold maximal mixing and twofold maximal z/^ — Vr mixing are 
consistent with the combined data then, it is clear that in threefold maximal mixing 
the best-fit value of Am^ is extremely close to the published CHOOZ limit |]T| on Am? 
computed for the case of twofold maximal — z/g mixing. This near-coincidence of 
the best-fit value with the published upper-limit strongly suggests that if threefold 
maximal mixing is valid, either of the two existing long-baseline reactor experiments 
CHOOZ |T| or PALO-VERDE |2^ could, with continued running, increased detector 
mass and especially with increased baseline L ~ 5 — 10 km, still be the first exper- 
iments to observe 'man-made' neutrino oscillations, and to establish the existence of 
Vf. mixing at the 'atmospheric scale' Am? ~ 10"'^ eV^. Otherwise the KAMLAND 
experiment |2^ should be decisive. 



7. Perspective and Future Prospects 

Naturally there are many other mixing schemes which can similarly describe the data. 



We mention explicitly here the Fritzsch-Xing hypothesis |30|, and in particular also 



the 'bi-maximal' hypothesis ||31[ and its generalisations |^3[ in which the is 
assumed to have Vei ^^ii content 0, 1/2, 1/2 perfectly replicating the phenomenology 
of twofold — mixing as far as the atmospheric experiments are concerned. In 
such schemes the zero or near-zero in the top right-hand corner of the mixing matrix 
guarantees no mixing at the atmospheric scale Am^ ~ 10^'^ eV^. 

Clearly certain 'bi-maximal' schemes can mimic 'tri-maximal' mixing (ie. threefold 
maximal mixing) extremely effectively. For example a bi-maximal scheme with the 
z/2 (tri-) maximally mixed (ie. having i/g, ^'/x, t-'t content 1/3, 1/3, 1/3), would have 
very similar phenomenology to threefold maximal mixing even for the solar data, with 
P(e — > e) = 5/9 in the gallium experiments and with the added possibility to exploit a 
large angle MSW solution with P(e — e) = 1/3 in the 'bathtub', in the higher energy 
experiments, for some particular Am'^ ~ 10^^ eV^. We emphasise however that in the 
threefold maximal mixing scenario suppression of Ve mixing at the atmospheric scale 
Am^ ~ 10~^ eV^, occurs naturally as a result of terrestrial matter effects without the 
insertion of an arbitrary zero in the mixing matrix as in the bi-maximal schemes. 

As regards 'bi-maximal' versus 'tri-maximal' mixing then, the crucial experimental 
question would seem to be whether or not there is appreciable vacuum v^, mixing at 
the atmospheric scale. This question might perhaps be answered in the future by 
the SUPER-K experiment itself, with the observation of a significant upward excess 
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of e-like events in the atmospheric data at muhi-GeV energies. The predicted U/D 
ratio for v^. in threefold maximal mixing, when matter effects are included, peaks as a 
function of neutrino energy at ~ 7 GeV, where {U / D)^ ^ 1.13. 

Together with reactor experiments, long-baseline accelerator experiments using 
high-energy z/^-beams will also probe z/g mixing. Fig. 8 shows the predicted energy- 
averaged appearance probability < P{fi ^ e) > for various mean neutrino energies in 
threefold maximal mixing for Am^ = 1.0 x 10^^ eV^. The curves are calculated with 
and without matter effects (solid and dotted curves respectively) and for visual clarity 
are shown averaged (50 : 50) over z//z/ beams. The suppression of z/g appearance as a 
result of matter effects is plainly visible at large L, but for all proposed experiments 
with L ^ 1000 km, matter effects can be largely neglected. 



8. Conclusion 



In this paper we have shown that the recent data from CHOOZ and SUPER-K 
are consistent with threefold maximal mixing (ie. 'tri-maximal' mixing) for a neutrino 
mass-squared difference Am^ ~ (0.98 ± g'23) ^ 10~^ eV^. With just a very few 
acknowledged exceptions, notably the LSND appearance result and possibly also 
the HOMESTAKE solar measurement Jl^ , the totality of previous data are also known 
P to be in good agreement with this hypothesis 0. 

Assuming tri-maximal mixing therefore, and taking a hierarchical spectrum of neu- 
trino masses similar to that of the quarks and leptons, we have for the heavy neutrino 
mass m3 ~ 30 ± 5 meV. This result is consistent with a 'see-saw' relation to the up- 
type quarks for a heavy (right-handed) Majorana mass ~ 1.0 x 10^^ GeV. The 
compton wavelenth of the heavy neutrino z/3 is A3 ~ 1/20 mm. Such a neutrino would 
make little contribution to cosmological dark matter. 
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3-Fold Max. (mat.) 
Am^ = 0.98 X 10-3 eV^ 


2-Fold Max. (i/^ - i/^) 
Am^ = 2.2 X 10-=^ eV^ 


Zenith-Angle 
Multi-GeV 


xVDOF = 6.5/7 
(CL = 48%) 


xVDOF = 2.0/7 
(CL = 96%) 


+ 

Zenith- Angle 
Sub-GeV 


xVDOF = 11.9/15 
(CL = 69%) 


xVDOF = 7.1/15 
(CL = 95%) 


+ 

R- Values 
Mult+Sub 


xVDOF = 13.2/17 
(CL = 72%) 


xVDOF = 7.6/17 
(CL = 97%) 


+ 

Upward- Muons 
Stop/Thru 


xVDOF = 13.6/18 
(CL = 75%) 


xVDOF = 7.7/18 
(CL = 98%) 


+ 

CHOOZ 
P(e ^ e) 


xVDOF = 25.4/28 
(CL = 61%) 


xVDOF = 18.7/28 
(CL = 91%) 



Table 1: Cumulative breakdown of x^ contributions and confidence levels for the 
various data-sets (see text) calculated for threefold maximal mixing with Am^ = 
0.98 X 10"^ eV^ and twofold maximal — mixing with Am^ = 2.2 x 10"^ eV^. The 
Am^ values quoted correspond to the overall x^-v[Amm\xm in each case summed over 
all the data-sets listed. Both threefold maximal mixing and twofold maximal i/^ — 
mixing give excellent overall fits to the data. 
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Figure Captions 



Figure 1. a) The combined KAMIOKA/SUPER-K multi-GeV zenith angle distribu- 
tion for a) /x-hke events and b) e-hke events. The sohd curve in each case is the 
threefold maximal mixing prediction for Am^ = 0.98 x 10^^ eV^, including terres- 
trial matter effects and compensation effects. The twofold maximal z/^ — mixing 
predictions (dashed curve) and the threefold vacuum predictions (dotted curve) are 
also shown. Either threefold maximal mixing with matter effects or twofold maximal 
— mixing are consistent with the multi-GeV zenith-angle data. 

Figure 2. The combined KAMIOKA/SUPER-K data on the up/down ratios for e- 
like and //-like events with | cos^| > 0.2 (data-points and shaded bands). The various 
theoretical predictions are plotted versus iSrv? . Threefold maximal mixing with matter 
effects (solid curves) approaches two-fold maximal — mixing (dashed curves) as 
Am^ 0, and over-compensation of the i/g-rate (expected in vacuum, dotted curve) 
is suppressed, consistent with the data. 

Figure 3. The SUPER-K sub-GeV zenith angle distributions (lepton momentum 
p > 400 MeV) for a) //-hke and b) e-hke events. The curves correspond to threefold 
maximal mixing with terrestrial matter effects and compensation effects for Am^ = 
0.98 X 10"^ eV^. The dashed curve shows the expected energy-independent 'matter 
oscillations' of amplitude ±1/6 and density-dependent wavelength, which are obscured 
by angular smearing (solid curve) in the water-Cherenkov experiments. 

Figure 4. The atmospheric neutrino ratio R = {n/e)oATA/{f^/G)MC as measured by 
SUPER-K (data-points and shaded bands). Within the statistical and systematic 
errors shown the data for R are consistent with Am^ ~ 10^^ eV^. Threefold maximal 
mixing including terrestrial matter effects (solid curves) predicts less change in < i? > 
from sub-GeV to multi-GeV energies (open and filled data points respectively) than 
twofold maximal i/^ — i/j- mixing (dashed curves). 
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Figure 5. The double ratio S of stopping to through muons in data and monte-carlo 
in SUPER-K (data point and shaded band). Threefold maximal mixing with matter 
effects included (solid curve) predicts a somewhat lower minimal suppression, but is 
otherwise very similar to the twofold maximal v^j^—Vj. mixing prediction (dashed curve), 
with the shift of a factor 2/3 on the Am^-scale, see Eq. (14). The data for S point to 
Am^ ~ 10~^ eV^ in both twofold and threefold mixing. 

Figure 6. The data from the CHOOZ reactor experiment on the Ve survival probability 
P(e — > e) plotted as a function of neutrino energy E = 1.8 + where Ee is the mea- 
sured positron energy. The solid curve shows the threefold maximal mixing prediction 
for Am^ = 0.98 x 10~^ eV^. With continued running, increased detector mass and 
particularly with increased baseline L ~ 5 — 10 km, significant Pg disappearance might 
well be established in reactor experiments. 

Figure 7. The overall in our combined fit, plotted versus Am^. In twofold maximal 
T^^i — mixing (dashed curve) there is a broad minimum extending over the range 
Am^ ~ 7 X 10"'' — 4 X 10"''^ eV^ with a best-fit Am^ ~ 2.2 x 10"^ eV^ corresponding 
to an excellent fit (x^/DOF = 18.7/28, CL = 91%). In threefold maximal mixing 
(solid curve) Am^ is relatively precisely determined: Am^ ~ (0.98 ±0^23) x 10""^ eV^, 
but still corresponds to a very good fit (x^/DOF = 25.4/28, CL = 61%). 

Figure 8. Suppression of i/g mixing by terrestrial matter effects. The predicted 
appearance probability as a function of the propagation length L through the Earth, 
in threefold maximal mixing with Am^ = 10~^ eV^. The dotted curve is for vacuum 
mixing only and the solid curve includes matter effects (for visual clarity the curves 
plotted refer to a 50 : 50 average over u/u beams). Matter oscillations, with reduced 
wavelength {^/2GNf,)~^ , are evident at large L, but proposed long-baseline accelerator 
experiments with L ^ 1000 km are largely unaffected. 
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